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SUMMARY 

A method baaed on the concept of Station Functions is pre- 
sented for calculating the modes and the frequencies of nonuniform 
cantilever beams vibrating in torsion, bending, and coupled 
bending-torsion motion. The method combines some of the advan- 
tages of the Rayleigh -Ritz and Stodola methods, in that a con- 
tinuous loading function for the beam is used, with the advantages 
of the influence-coefficient method,- in that the continuous load- 
ing function is obtained in terms of the displacements at a finite 
number of stations along the beam. 

The Station Functions were derived for a number of stations 
ranging from one to eight. The deflections were obtained in terms 
of the physical properties of the beam and Station Numbers, which 
are general in nature and which have been tabulated for easy refer- 
ence. Examples were worked out in detail; comparisons were made 
with exact theoretical results. For a -uniform cantilever beam, for 
n .stations along the beam, the first n-1 modes and frequencies were 
in excellent agreement with the theoretically exact values. It was 
shown that the effect of coupling between bending and torsion is to 
reduce the first natural frequency to a value below that which it would 
have if there were no coupling. 


INTRODUCTION 

The failure of turbine and compressor blades due to vibrations 
has led to an increased interest in the study of the vibrations of 
these blades and in the determination of the natural modes and fre- 
quencies. In such theoretical studies, it is usually assumed that 
the compressor or turbine blade acts as a cantilever beam. The 
calculation of the uncoupled modes of arbitrarily shaped cantilever 
beams has been extensively investigated (references 1 to 4), but 
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little work has as yet been done on calculating the coupled modes 
of such "beams. If the geometry of the beam is such that coupling 
exists, the coupled modes are the actual vibrational modes that 
must be calculated. 

Four general methods are currently in use for calculating 
uncoupled modes and frequencies of nonuniform beams . These methods 
are the Rayleigh-Bitz or energy method (reference l), the Stodola 
method (references 5 and 6), the influence-coefficient method (refer- 
ences 4 and 7), and the integral -equation method (references 8 
and 9). For each of these methods, computational work can usually 
be carried out in several ways. For example, by the use of influ- 
ence coefficients the modes and frequencies can be determined by 
Mykelstad’s iteration procedure (reference 7) or by matrix 
methods (reference 4). 

Any one of these methods can be extended to the calculation 
of coupled bend ing-tors ion modes. The Rayleigh-Bitz method usually 
requires that the uncoupled modes be determined before the coupled 
modes can be computed. In applying either the Rayleigh -Bitz or 
the Stodola method, great difficulty is encountered in accurately 
determining the higher modes, because the lower modes must first 
be "swept out” by the use of exact orthogonality conditions (refer- 
ence 10); the process will otherwise always converge back to the 
lowest mode. The same difficulties are encountered in the integral- 
equation method. 

The influence-coefficient method avoids these difficulties by 
reducing the problem to one having a finite n umb er of degrees of 
freedom. The beam is divided into n intervals and a concentrated 
loading is assumed at the center of gravity of each interval. The 
solution of the resultant determinantal equation gives the first 
n modes. The accuracy of the higher modes is, however, very poor. 
Only the first third of the modes and the first half of the fre- 
quencies are obtained within the usual engineering accuracy. 

Carrying along so many useless modes greatly increases the labor 
involved. 

A straightforward, accurate method for determining the coupled 
bending-torsion modes and the frequencies of nonuniform cantilever 
beams, together with applications of this method, was developed at 
the NACA Lewis laboratory and is presented herein. This method is 
based on the use of Station Functions as first discussed in refer- 
ence 11. Incorporated in the method are the advantages of the 
continuous -function deflections of the Rayleigh-Bitz and Stodola 
methods together with the advantages of the finite number of 
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degrees of freedom of the influence-coefficient method. When the 
method is applied to a uniform beam, the first n-1 roots of the 
resultant determinantal equation are amply accurate for engineer- 
ing purposes. 

The final determinantal equation is solved herein by matrix- 
iteration methods. Any other convenient method may, however, be 
used and no knowledge of matrix algebra is needed to carry out the 
calculations by the matrix method. The work can be done by an 
inexperienced computer, as the only operations necessary for 
determining each mode are cumulative multiplication and division. 

In addition, for the case in which the coupling coefficient remains 
constant along the beam, a simple quadratic formula and a series of 
curves are presented for determining the first coupled mode in 
terms of the uncoupled modes . Examples are developed in detail 
and comparisons with exact theoretical results are included. 


THEORY 

In the usual inf luence -coefficient methods for solving dynam- 
ical problems, a continuous body having an infinite number of 
degrees of freedom is replaced by a body having a finite number of 
degrees of freedom. Two principal assumptions are then made, which 
introduce inaccuracies into the solutions, particularly in the 
higher modes: (1) The resultant of the inertia loads of all the 

infinitesimal masses in a finite interval passes through the center 
of gravity of that interval; and (2) a concentrated load that is 
the resultant of a distributed load produces the same deflection 
as the distributed load. An attempt has-been made to reduce the 
error due to the second of these assumptions by the use .of weight- 
ing matrices (reference 12). Although the accuracy is thereby 
increased, the effect of the first assumption is still great enough 
to introduce serious errors (reference 11). 

In order to eliminate these assumptions, Rauscher (refer- 
ence 11) introduced the concept of Station Functions. Instead of 
assuming the inertia loads to be concentrated at the centers of 
gravity of the intervals, the inertia loads and, consequently, the 
deflections are assumed to be continuous functions along the beam. 
The values of these continuous deflection functions at the refer- 
ence stations must equal the deflections of the reference stations. 
The loading on the beam is therefore a continuous function of the 
deflections of the reference stations. Inasmuch as the deflections 
of the reference stations can be computed from the loading on the 
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■beam, which in turn is available from the deflections, the deflec- 
tions are therefore obtained as functions of themselves. This pro- 
cedure gives n homogenous equations in the n deflections of the 
reference stations. The resultant determinantal equation has 
n roots' for the frequency; it will be shown that at least n-1 of 
these roots are sufficiently accurate for engineering purposes if 
the deflection functions are properly chosen. (For coupled bending- 
torsion vibrations, 2n homogeneous equations and 2n roots are 
obtained for n stations . ) 

The deflection functions used must satisfy the boundary con- 
ditions of the problem and also the condition that, at any refer- 
ence station, the value of the function must equal the deflection 
of the reference station. Although it is always possible to find 
directly a single function that will satisfy these conditions, it 
is more convenient to obtain different component functions at each 
station and to add all these component functions together to give 
the complete deflection function. Rauscher (reference 11) calls 
these component deflection functions Station Functions. For 
example, the complete torsional deflection function for the beam 
will have the following form: 

0 U) = £ fiU) . 

j=l 3 3 

where 

z dimensionless distance along beam 

0(z) torsional deflection at distance z from root 

torsional deflection at the station • 

f j(z) Station Function in torsion associated with the 

station (All symbols are defined in appendix A.) 

Each Station Function must satisfy the boundary conditions of 
the problem and the following additional conditions: (1) At the 

reference station with which it is associated, the Station Function 
equals the deflection of that reference station; and (2) at all 
other reference stations, the Station Function equals zero. The 
sum of all these Station Functions will then give the complete 
deflection function for the beam. The Station Functions and cor- 
responding loading functions are derived in appendix B for tor- 
sional vibrations, bending vibrations, and coupled bending-torsion 
vibrations of an arbitrary cantilever beam. 
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Torsional vibrations . - It is shown in appendix B that the 
torsional deflections of the reference stations for a team divided 
into n intervals of length 6, as shown in figure 1, are given 
"by the following system of equations : 


0> 2 5 2 ija 

°o6i 


ij e i 


where 


k 

i and J = 1,2, . . . n 


ij -Z. — | A k N Jk - ( k ’ 1 ) h M Jk + ZZ I r M jr 
k=l ^k L r=k+l 


CO frequency of vibration 

& length of interval 


( 1 ) 


( 2 ) 


1(3 mass moment of inertia per unit length about elastic axis at 
root section 

Ifc ratio of average mass moment of inertia per unit length of 

k th interval to mass moment of inertia per unit length at 
root section 


Cq torsional stiffness of the root section 

C k ratio of average torsional stiffness of k^* 1 interval to 
torsional stiffness at root section 


The Station Numbers 
integers k, J, and n 


Njjj- and Mjjf are functions only of the 
and are defined as 


N Jk S 


yj 


Pk 

k-1 

Pk 


"\ 


Zf j(z)dz 


(3) 


Hjk* 


yj 


k-1 


f j(z)dz 


J 


where f i(z) represents the Station Functions derived in appendix B 
and is given by 
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fj(z) « ay z + ay z 2 + . . . +• a (n+1)J z( n+1 ) (4) 


The coefficients ay are determined in appendix B by satis- 
fying the conditions on the Station Functions . The integrals in 
equations (3) are thus seen to he integrals of simple polynomials 
and the limits of integration are integers. The Station Num- 
bers Njk and Mjk are therefore rational numbers, functions only 
of the integers n, k, and J . These numbers have been evaluated and 
and. are listed in tables I to VIII. 


If the physical properties of the beam under consideration sire 
known for each of the n intervals, and 1^ will be known. 

The Station Numbers Nj^ and can be obtained from tables I 

to VIII. From equation (2), ay can then be easily calculated. 


Equation (l) actually represents n homogeneous equations in 
the n unknown deflections 0y With i_ = x, these equa- 
tions can be written as follows : w 

(an-x) e 1 + a 12 e 2 + a 13 e 3 + . . . + a ln e n = 0 ^ 

a 21 0 1 + ( a 22‘-^) e 2 + ^Z 9 ^ + • • • + a 2n e n “ 0 
a 31 9 1 + .°32 °2 + ^ZZ~^ e 3 + • • • + a 3n e n ■ 0 \ (5) 

“hi 0 ! + °n2 0 2 + “nS 0 3 + • • • + (“jan"*) e n ■ 0 j 


For a nontrivial solution, the determinant of the coefficients must 
vanish and the characteristic equation becomes 


a ll-X 

a 12 

a 13 • • 

• a ln 



a 21 

a 22“^ 

a 23 * • 

• a 2n 



a,, 

31 

a 32 

a 33 -X . . 

* a 3n 

« 0 

( 6 ) 

• • • 

“nl 

• * • 

a n2 

“113 * * 

• « • 

• “on-* 
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or 

XI - pij] =0 (6a) 

where I is the Identity matrix, and [aijj is the dynamical 
matrix. 

Equation (6) can he solved for the n values of X by any 
method available. The method used herein was to obtain the values 
of X as the latent roots of the matrix [a^j] , which is actually 
the dynamical matrix for the problem. The mode shapes are obtained 
at the same time. 

Bending vibrations . - The bending deflections for the beam 
shown in figure 1 are given by the following system of equations 
(see appendix B): 


04 ®0 

y i 8 BX 


^ V a 

B 0 4 =* 


where 






i and J = 1,2, . . . n 

m Q mass per unit length of beam at root section 

ratio of average mass per unit length of k^* 1 interval to 
mass per unit length at root section 

Bq bending stiffness at root section 

Bfc ratio of average bending stiffness of k^ interval to bend- 
ing stiffness at root section 
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The Station Numbers M ' j k , N' j k , 
functions only of the integers k, j, 

r-k 


P'jk, and Q» Jk are 
and n and are defined by 



(k-1) 



gj(z)dz 




^k-1 


z gj(z)dz 


J 


The Station Functions gj(z) are derived in appendix B and are 
given by 

g J (z) ’ b 2J ** + b 3J * 3 * b 4J 2 * + • • • + b (n+3)4 2<n+3) 


( 10 ) 

The integrals in equations (9) are thus seen to be integrals of 
simple polynomials. The Station Numbers M'j k> N' j k , P'jk* 
and Q' jk are rational numbers, functions only of the integers j, 
k, and n. These numbers have been evaluated and are listed in 
tables I to VIII. 

If the physical properties of the beam are known for each of 
the n intervals, m k and Bk will be known. The Station Num- 
bers N'jk* -P'jk* and Q'jk 611,6 obtained from tables I 

to VIII; 3jj can then easily be calculated by equation (8). 


The determinantal equation is 



